ADA  1 10440 


i  MfM 

&  '*f4  jfe  y* 

rlNi  ii'  k*®® 


NCL  Memorandum  Report  4727 


Surface  Effects  Due  to  Subsurface  Processes: 

A  Surrey 


D.  T.  ClIEN 

Space  Serving  Applications  Branch 
Aerospace  Systems  Division 


January*!  5,  1982 


DTfC 

ELECTEI 
FEB  4  W82 


mSW 


NAVM,  *FSF.\»r  ?!  !.  VBORATO*  Y 
FVash^ne^rt.  DC 


StCUAlTv  CLASSIFICATION  OF  This  FAOt  rWiAn  Omlm  En m.aj 


REPORT  DOCUMENTATION  PAGE 


9>/;e  READ  INSTRUCTIONS 

BEFORE  COMPLETING  FORM 
2.  OOVT  ACCESSION  NO.  1  AEClAlENT'S  CATALOG  NUMICA 


NRL  Memorandum  Report  4727 


P-iAilo 


4.  TITL t  JuAlllU; 

SURFACE  EFFECTS  DUE  TO  SUBSURFACE  PROCESSES: 
A  SURVEY 


S.  YYAE  OF  AEPOAT  4  FIMOO  COVERED 

Interim  report  on  e  continuing  NRL 
proble.a. 

(.  PERFORMING  oao.  acaoat  number 


IT.  AUTMOAf*; 


D.  T.  Chen 


(.  PERFORMING  OAGANIIATION  NAME  AnO  ADDRESS 

Naval  Research  Laboratory 
Washington,  DC  20375 

M.  CONTROLLING  OFFICE  NAME  ANO  ADDRESS 

Office  of  Naval  Research 
Arlington,  VA  22217 


I.  CONTRACT  OA  OAANT  nuMBERC.; 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AACA  A  WOAK  UNIT  NUMBERS 

61153N;  RR032-08|-41;  79-0930-0  A 
61153N;  RR032-04^02;  79-C001-X 

12.  ACAOAT  DATE 

January  15, 1982 

I).  NUMBCK  OF  RAGES 

44 


4.  MONiTOAInG  AGENCY  name  4  ADDRESS//!  dlliitmnt  Horn  Conlnlllnt  Olttc •)  IS.  SCCUAITV  CLASS.  Co/  Mi/a  nfMj 

UNCLASSIFIED _ 

IS  A.  OCCl  assi  fi  Cation/  lownoraOinq 
schedule 


IS-  OISTAISUTIOH  STATEMENT  (ol  thl,  rt.por r 


Approved  for  public  release;  distribution  unlimited. 


[  17.  DISTRIBUTION  STATEMENT  (of  thb  mbatrmet  mnfrmd  In  Block  20,  It  dttloront  trom  Report) 


ha.  supplementary  NOTES 


1 19.  KEY  WORDS  fContlnuo  on  tormram  wldo  1 1  nm cooaory  Identify  by  Jtock  numb*) 


Weak  interaction# 

Strong  interactions 
Bathymetry 

Subsurface  meso scale  features  '  S' 
Internal  waves 


Subsurfaio  Ekman  layer 
Breaking  waves 
Significant  slope 

f  t'  r  V  '  . 
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The  purpose  of  this  document  is  to  provide  e  survey  oythe  oceanic  subsurface  dynamic  pro¬ 
cesses  which  have  signatures /or  affactjyon  the  ocean  surface  These  surface  effects  can  appear  in  the 
form  of  changes  in  both  surface  current  distributions  and  local  nave  structures  (or  profiles).  Based 
upon  their  physical  scales,  these  effects  are  'OBtegepil'*- classified  into  phenomena  p  reduced  by  mech¬ 
anisms  associated  with  weak  and  strong  interactions.  Two  types  of  phenomena  are  distinct  among 
surface  wave  weak  interactions.'  The  -first  typ^broa  d  band  phinomeri*^  includes  wa^wave  energy 
transfer  and  general  wave  train  instability.  TS*  second  type^  narrow  band  phenomena'^  includes 
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B*njaminFe  ^instability,  neumnet,  and  envelop*  solitons.  All  these  phenomena  riser*  the/physical 
characteristic^®* the  following  (1)  they  may  be  manifested  only  by  the  long  dominant  waves  in  a 


wave  field  and  yet  the  behave  of  riiort  waves  is  determined  by  the  strong  interactions  with  the  long 
wave*  rather  than  by  these  ^roceeses;>(2)  they  are  evolutionary  phenomena  with  a  time  scale  of 
rt*h)^j  tiiws  the  wave  period  where  a'  is  the  wave  amplitude  and  k  is  the  wave  number  of  the  dom- 
TSSai  wavs.  frhtwaaa^jw  strong  interactions  are  those  phenomena  whose  time  scale  is  of  the  order  of 
the  wav*  period  and  space  scale  is  of  the  order  of  the  wavelength.  JFh^de tailed  wav*  structures  (or 
profiles)  are  the  prime  sources  of  information.  Strong  interactions  include  the  strong  Longuet-Higgina 
instability,  wave-current  interactions,  long  wave-short  wave  interactions,  and  the  processes  of  parasitic 
capillary  formation  and  microscale  breaking  induced  by  surface  wind  drift  jADue  to  the  rapid  surface 
wave  responses  to  local  oceanic  disturbances  accounted  by  the  strong  intemtiona,  their  phenomena 
can  offer  at  least  qualitatively  aathe  indicator  of,  for  example,  the  exiatenosW  internal  wave 
patterns.  t  \ 

Surface  effects  due  to  bathymetry,  subsurface  mesoacale  features,  internal  waves,  subsurface 
Ekman  layer  and  breaking  wave,  and  their  relationship  to  the  subsurface  processes,  then,  are  discussed 
within  the  context  set  forth  for  weak  and  strong  interactions.  Finally,  in  order  to  highlight  the 
drain-. ic  recent  scientific  and  technical  advances  made  for  the  space  age,  a  new  approach  to  simulate 
a  wide  variety  of  both  surface  and  subsurface  dynamic  phenomena  is  briefly  mentioned  for  evaluating 
the  influence  of  breaking  waves  on  the  dynamics  of  the  upper  ocean  by  using  a  single  measurable 
parameter,  significant  slope,  from  a  spacebome  altimeter. 
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SURFACE  EFFECTS  DUE  TO  SUBSURFACE  PROCESSES:  A  SURVEY 


1 .  INTRODUCTION 

The  purpose  of  this  document  Is  to  provide  a  survey  of  the  oceanic 
subsurface  dynamic  processes  which  have  signatures,  or  effects,  on  the  ocean 
surface.  These  surface  effect*  can  appear  In  the  form  of  changes  In  both 
surface  cur.-ent  distributions  and  local  wave  structures  (profiles).  Although 
the  microscale  breaking  of  the  surface  wave  disrupts  the  diffusive  surface 
layer,  the  thermal  surface  effects  are  not  Included  because  of  the  uncertainties 
In  the  mechanisms  Involved  In  the  heat  transfer  across  the  air-sea  Interface. 
Interested  readers  should  refer  to  Phillips  (1974)  on  the  subject  of  thermal 
surface  effects  or  of  slicks.  Nevertheless,  the  thermal  surface  effects  are 
treated  here  only  as  the  surface  temperature  boundary  conditions  in  the  context 
of  oceanic  transport  of  heat. 

Oceanic  subsurface  dynamic  processes,  in  this  case,  can  br  those  shown  In 
the  Internal -wave-modi fled  current  and  wave  field,  the  mesoscale  subsurface 
eddies  and  fronts,  and  the  bottom- topography-modi  fled  current  field.  Except 
for  the  Internal  wave  Induced  field,  which  has  an  additional  direct  wavy 
surface  pattern,  all  the  subsurface  dynamic  processes  have  surface  current 
effects.  These  surface  current  effects  present  themselves  as  variations  in 
amplitudes  and  directions  over  the  underlying  current  field.  The  magnitude  of 
these  variations  can  be  of  the  order  of  cm/sec,  i.e.,  small,  but  their  effects 
can  be  manifested  through  wave-current  and  wave-wave  interactions.  Although 
it  is  not  the  intent  of  this  document  to  identify  the  remote  sensing  techniques 
by  which  the  variations  in  current  patterns  and/or  the  modulation  of  surface 
waves  can  be  observed  (it  should  be,  however,  the  subject  of  the  next  document), 
it  is  easier  to  observe  remotely  the  modulation  of  surface  wavec.  Therefore, 
this  document,  more  or  less,  concentrates  on  the  response  of  the  wave  field  to 
a  moving  current  pattern  or  stationary  current  variations  generated  by  subsurface 
dynamic  processes. 

2.  SURFACE  WAVE  EFFECTS:  THE  NATURE  OF  WEAK  AND  STRONG  INTERACTIONS 

It  has  been  known  that,  surface  waves  of  sufficiently  long  wavelength  and 
sufficiently  small  slope  can  propagate  for  long  distances  without 
significant  change  of  their  wave  form.  This  phenomenon,  basically,  asserts 
the  utilization  of  the  principle  of  superposition,  whose  validity  requires 
linear  relationship  without  mutual  interference.  In  waves  with  sufficiently 
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small  slope.  However,  these  waves  have  also  been  observed  to  be  dispersive  In 
amplitude  as  well  as  In  phase  after,  say,  a  time  of  100  wave  periods  or  a 
distance  of  100  wavelengths.  It  is  obvious  that  there  are  small  regular 
perturbations  not  only  In  amplitude  but  also  in  phase  to  the  linear  solution 
of  the  governing  equations.  The  outstanding  response  of  the  waves  to  these 
small  perturbations  is  slow  and  small.  The  effect,  which  may  be  cumulatively 
large,  cannot  be  appreciated  significantly  In  local  space  and  time  scales.  In 
fact,  locally  the  waves  are  almost  like  linear  waves,  but  as  a  result  of  the 
small  regular  perturbations,  they  change  very  slowly  In  space  or  time  or  both. 
These  changes  are  classified  as  the  weak  Interactions  which  Include  wave-wave 
Interactions  (Phillips,  1977),  Benjamln-Felr  instabilities  (Benjamin  and  Felr, 
1967),  the  existence  of  envelope  solltons  (Yuen  and  Lake,  1975),  etc.  It  Is 
the  nature  of  the  changes  which  Is  of  the  third  order  and  is  of  Interest,  not 
the  Individual  character  of  each  wave.  Except  for  the  bottom-bathymetry 
induced  wave-wave  Interactions  (Herterich  and  Hasselmann,  1980),  It  should  be 
noted,  however,  that  all  the  weak  interactions  are  surface  wave  induced  phenomena. 

In  contrast  are  the  strong  Interactions  which  can  be  interaction  of  a 
wave  with  itself,  another  wave  or  a  current.  The  changes  are  local  occurrences 
but  widespread  and  attributed.  The  time  scale  of  these  changes  is  comparable 
with  or  shorter  than  the  wave  period.  Therefore,  the  individual  structures  (or 
profiles)  of  the  waves  are  important  quantities  to  be  measured.  Strong  inter¬ 
actions  include  the  phenomenon  such  as  wave  breaking  with  the  formation  of 
whitecaps.  Wave  breaking  certainly  Is  a  non-linear  phenomenon  and  cannot  be 
treated  mathematically  as  a  small  perturbation  to  the  linear  solution. 

The  difference  between  the  strong  and  the  weak  Interaction  phenomena  is 
distinctly  fundamental:  weak  Interaction  phenomena  take  many  wave  periods  or 
many  wavelengths  to  become  signficant  although  their  cumulative  effects  may  be 
large  whereas  strong  interaction  phenomena  occur  rapidly  and  locally  although 
their  net  effects  may  be  small.  Therefore,  wave  property  changes  locally  should 
be  accounted  for  by  the  strong  interactions,  and  wave  property  changes  in  a 
gradual  manner  over  a  long  distance  or  long  duration  of  time  should  be  accounted 
for  by  the  weak  interactions.  Phenomena  produced  by  the  weak  Interactions  are 
considered  as  noise  when  one  is  seeking  to  account  for  phenomena  produced  by 
the  strong  Interactions,  and  vice  versa.  Although  this  document  is  primarily 
Interested  in  the  wave-current  interaction— a  strong  interaction— all  aspects 
of  the  relevant  weak  and  strong  interactions  should  be  explored. 
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3.  WEAK  INTERACTIONS 


The  essential  approximation  in  the  study  of  the  weak  interaction  is  that 
the  wave  shape  ak  is  small  where  a  is  wave  amplitude,  k  ■  2ir/L  is  wave  number, 
«md  L  is  wavelength.  Let  ;  be  the  wave  surface  displacement  at  position  x 
and  f-*..»e  t  with  respect  to  the  mean  sea  surface  along  the  vertical  axis  z, 
positive  upward.  The  nonlinear  boundary  conditions,  specified  at  the  free 
surface  z  ■  c(x,  t),  can  be  referred,  by  means  of  a  Taylor  series  expansion, 
to  the  mean  sea  level  z  *  0  and  the  resulting  terns  can  be  collected  In  the 
order  of  ak.  This  is  the  Stokes'  expansion  and  its  Important  role  In  the 
standard  procedure  for  the  solution  of  a  set  of  nonlinear  governing  equations 
for  non-zero  but  small  ak  is  well  known.  Restriction,  commented  by  Phillips 
(1979),  on  the  small  amplitudes  of  the  longest  waves  imposed  by  Stokes' 
expansion  about  the  mean  sea  level  z  *  0  can  be  alleviated  by  expanding  on  the 
longest  wave  surface  with  the  Lagranglan  reference  frame  which  moves  with  this 
longest  wave  surface  (Liu  and  Benney,  1980). 

Employing  the  techniques  of  Stokes'  expansion,  theoretical  studies 
(Phillips,  1977)  indicate  that  energy  exchange  and  amplitude  modulation  of 
weak  interactions  appear  at  the  third  order,  not  at  the  first  and  the  second 
orders,  and  only  among  sets  of  wave  components  whose  wave-number  vectors  form 
a  closed  quadrilateral,  which  satisfies, 

+  ^2  3  ^3  +  k4  (3.1) 

and  whose  wave  frequencies  satisfy  the  resonant  conditions, 

al  +  °2  3  °3  +  °4  (3.2) 


where 


Oj  3  glfcjl.  i  3  1.  2,  3,  4  (3.3) 

g  is  the  gravitational  acceleration  and  |tj|  is  the  magnitude  of  the  wave- 
number  vectors  It ^ .  Because  theoretical  results  predict  that  terms  describing 
energy  exchanges  and  amplitude  modulations  are  of  the  third  order,  the  time  or 
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.p 

spatial  scales  for  weak  Interactions  are  (ak)  times  the  wave  period  or 
wavelength,  respectively. 

The  governing  equations  for  weak  Interactions  are  found  to  be 


• 

al  al 

al911  alg12  algl3  alg14 

★  ' 
a2  a3  a4 

o2  a2 

a2g21  a2g22  a2g23  a2g24 

.  J 

1a2*Z 

al  a3  a4 

< 

• 

►  < 

* 

}  +  1h  \ 

°3  *3 

a3g31  a3g32  a3g33  a3g34 

1c3*3 

a1  a2  a4 

V°4S4 

,a4g41  a4g42  a4943  a4g44 

ia4a4 

al  a2  a3; 

where  denotes  time  derivative  (3/3t),  designates  complex  conjugate,  1 

Is  the  Imaginary  number,  a^'s  are  wave  amplitudes,  g^  and  h  are  real  functions 

of  the  wave-number  vector  configuration,  represented  by  Equation  (3.1).  If 

wave  amplitude  varies  spatially  as  well  as  in  time,  a.  Is  replaced  by 

( 3/3t  +  c  »v)a<  where  c„  Is  the  group  velocity.  The  term  g. f  Is  the  square 
9  J  9 

matrix  In  Equation  (3.4)  and  describes  Interactions  of  wave  k^  to  wave  k^ 
per  unit  wave  energy,  and  gives  changes  in  phase  speeds.  These  changes  of 
phase  speeds  of  Infinitesimal  waves  associated  with  finite  wave  amplitude  are 
called  amplitude  dispersion  effect.  The  second  column  matrix  at  the  right 
hand  side  specifies  the  energy  exchange. 


The  phenomena  of  weak  interactions  Involve  the  interrelationship  of  the 
two  effects  -  amplitude  dispersion  and  energy  exchange.  Benney  (1962)  has 
shown  that  the  total  wave  energy  and  momentum  densities  of  each  tetrad  of 
waves  are  conserved.  The  general  solution  of  Equation  (3.4)  can  be  written  In 
terms  of  elliptic  functions  (Bretherton,  1964). 


3.1  Physical  Phenomena  of  Weak  Interactions 

3.1.1  Wave-Wave  Interaction 

Replacing  and  o^  in  Equations  (3.1)  and  (3.2)  by  t3  and  a^,  respec¬ 
tively,  Equations  (3.1)  and  (3.2)  reduce  to 

1 1  +  ]T2  -  21^3  (3.5) 
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and 


a-|  +  c2  3  2ct^  .  (3.6) 

The  existence  of  this  physical  phenomenon  of  energy  transfer  among  different 
wave  modes  was  confirmed  separately  by  Longuet-HIgglns  and  Smith  (1966)  and 
McGoldrlck  et  al.  (1966).  Because  the  experiment  performed  by  the  former  has 
considerably  higher  wave  slopes,  the  amplitude  dispersion  effects  were  signi¬ 
ficant  which  Implies  simple  shifts  In  the  resonance  band  to  the  smaller  wave- 
number  vector. 

3.1.2  The  Instability  of  a  Wave  Train 

When  there  Is  a  slight  mismatch  from  resonance  such  that,  from  Equation 
(3.6), 


al  +  a2”  ^a3  =  ^  *  (3.7) 

Then,  the  right  hand  side  of  Equation  (3.4)  contains  an  additional  factor 
e"^  (Longuet-HIgglns  and  Smith,  1966).  For  components  or  ic2  Equation 


(3.4)  would  be  reduced  to 

|  h2Oia2(A3A3)2o1  (3.8) 

“2“  ^  ^°i02^3^3^a2  (3.9) 

A3s  I933A3A3  (3.10) 

where 

a1  3  a]  exp(-1g13A3A*t)  (3.11) 

Og  3  a2  expHg^A-jAjt)  (3.12) 


and 


1 


as  wave  component  l<3,  represented  by, 

A3  *  «3('  +  0 

Is  perturbed  slightly  by  waye  component  or  ^  under  the  assumption  that  the 
mismatch  £  In  the  resonant  condition.  Equation  (3.7),  Is  compensated  by  the 
Influence  of  amplitude  dispersion,  l.e.. 
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-(9i3  +  923)A3A3 


(3.13) 


Solutions  for  and  from  Equations  (3.8)  and  (3.9)  are  unstable. 
Furthermore,  If  the  mismatch  £  cannot  compensate  the  amplitude  dispersion, 
as  shown  in  Equation  (3.13),  the  stability  analysis  Is  valid  only  for 
|a-j|  and  jagl  «  | [ •  The  detailed  description  of  the  Instability  of  a 
Wave  Train  should  be  found  in  Benjamin  and  Feir  (1967)  and  Yuen  (1977). 

3.2  Narrow  Band  Spectrum  Approximation 

For  a  wave  group  or  a  slowly  varying  wave  train  whose  surface  displace 

ment 


C(x,t)  =  exp[i(tQ*x  -aot)]Ja()t)exp{if!c*x  -  (o(|ic|  -ao)tj}dl<  (3.14) 

where  the  only  significant  contribution  come  from  a  small  range  of  wave- 
numbers  and  frequencies  aoout  t  and  oQ,  respectively,  Equation  (3.4)  can  be 
considerably  simplified  and  reduced  to  the  nonlinear  Schrodlnger  equation, 


I'?oIZooIaiHa  * 0 


(3.15) 


where  A  is  the  local  amplitude, 

A  a/a0<)  exp{i[£*x  -(o(|ltl)  -<70Jt])d£  (3.16) 

from  Equation  (3.14)  and  the  subscripts  t,  xx,  and  yy  represent  partial 


^1 
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derivatives.  This  equation  was  derived  by  Benney  ar.d  Newell  (1967).  The 
physical  phenomena  modelled  by  Equation  (3.15)  are  still  fundamentally 
those  of  weak  Interaction,  amplitude  dispersion  and  inter-mode  energy 
transfer. 

The  physical  dynamic  processes  of  Benjamin-Felr  Instability,  recurrence, 
and  envelope  solltons  can  all  be  classified  as  the  narrow  band  weak  Inter¬ 
action  phenomena.  Benjamin-Felr  Instability  is  a  special  case  of  the 
Instability  of  a  wave  train  (Benjamin  and  Felr,  1967).  Recurrence  Is 
essentially  a  phenomenon  predicted  by  numerical  integration  of  Schrfidlnger's 
equation.  Equation  (3.15).  Yuen  and  Ferguson  (1978)  have  shown  that  the 
solution  to  Equation  (3.15)  with  essentially  its  initial  condition,  will 
return  after  a  time  interval  of  order  (a 1 1c | )-^a-1  in  the  numerical  inte¬ 
gration  scheme.  There  are  some  questions  raised  against  the  validity  of 
this  numerical  solution,  for  example,  Huang  (1981)  stated  that,  from 
experimental  evidence,  the  wave  pattern  may  be  back  to  the  original  form 
but  the  wavelength  has  certainly  increased.  More  research  work  is 
required  in  order  to  resolve  the  differences. 

Envelope  soli  tons  are  groups  of  surface  waves  whose  envelope  propa¬ 
gates  without  change  in  the  wave  form.  Phase  shift  may  happen  if  the 
envelope  soli  tons  run  through  other  envelope  soli  tons,  but  their  basic 
forms  still  remain  the  same.  Yuen  and  Lake  (1975)  among  others  have 
given  the  exact  solution  to  the  one-dimens i oral  nonlinear  Schrddinger 
equation  for  sufficiently  small  slopes  and  band-widths.  Their  solution 
Indicates  the  involvement  of  a  balance  between  amplitude  dispersion  and 
linear  dispersion  about  the  carrier  frequency.  The  problem  of  stability 
of  an  envelope  soliton  to  infinitesimal  directional  perturbations  has 
been  studied  by,  for  example,  Saffman  and  Yuen  (1978). 

4.  STRONG  INTERACTIONS 

The  strong  Interaction  phenomena  In  waves  are  those  which  are 
produced  In  a  time  scale  or  spatial  scale  comparable  with  or  shorter 
than  the  wave  period  or  the  wavelength,  respectively.  Among  these 
phenomena  are  local  wave  breaking  (Longuet-HIgglns,  1978a  and  1978b; 
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Longuet-HIggins  and  Cokelet,  1978),  the  Interaction  of  waves  with  non- 
uniform  current,  and  the  ii.ix /action  of  waves  with  another  long  wave. 

These  phenomena  are  highly  nonlinear  In  general  and  are  the  result  of 
Interactions  with  a  finite  amplitude  disturbance. 

Due  to  the  rapid  response  of  the  strong  interactions  in  waves,  the 
Influence  of  weak  interactions  can  almost  be  ignored.  Because  of  the 
instrumentational  difficulties  in  observing  and  mathematical  difficulties 
in  analyzing  the  strong  interaction  phenomena,  they  are  much  less  explored 
than  the  weak  interaction  phenomena. 

4.1  Physical  Phenomena  of  Strong  Interactions 

4.1.1  Wave-Current  Interactions 

As  the  waves  consist  of  a  train  or  group  of  local  wave-number  It 
and  frequency  a,  neglecting  the  weak  interactions,  the  kinematic  conservation 
equation 

||  +  1  (j  +  It  •  ft)  *  0  (4.1) 

should  be  satisfied  where  ft’(x.t)  is  the  current  system  encountered  and 
a  ~  (g|t|)^.  The  distribution  of  energy  density  E,  in  this  case, 
should  be  governed  by  the  conservation  of  wave  action,  f. 

It  <£■)  +  $•  [<ff  +  V§]*°  (4-2) 

where  cg  is  the  group  velocity.  Equation  (4.2)  can  be  rewritten  as 

ft  An(§)  +  $  •  (O'  +  cg)  -  0  (4.3) 

where  f^  represents  the  rate  of  change  moving  with  the  energy  packet 
which  moves  with  the  combined  current  and  group  velocities. 

One  interesting  situation  is  that  as  the  current  distribution  moves 
through  tne  ambient  fluid  field  such  as  when  an  internal  wave  pulse  or 
train  moves,  the  phenomenon  of  blockage  (Longuet-Higgins  and  Stewart, 
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1961)  may  result.  Let  ft  be  the  velocity  of  the  current  pattern  and  the 
distribution  of  current  be  $  =  ft(x  -  |ft|t),  then  the  modulations  in 
surface  wave  energy  density  E  is  governed  by 

E i c  ]  ( +  y|c|  -  |ft|)  =  const.  (4.4) 

where  c,  ft,  and  ft  are  all  collinear  in  direction;  the  local  surface  wave 
phase  velocity  c  is  found  from  Equation  (4.1).  For  internal  waves,  in 
most  cases,  jft|  «  |c]  or  |ft|.  The  intensity  of  the  surface  wave  modulation 
becomes  large  as  the  magnitude  of  group  speed  ^cl  approaches  that  of 
the  propagation  speed  |ft|  of  the  pattern.  In  situtations  where  surface 
and  internal  waves  are  oblique  with  each  other,  the  treatise  is  the  same 
as  for  the  collinear  situation  except  the  velocities  c,  ?  and  ft  are 
projected  in  the  same  direction  (Longuet-Higgins  and  Stewart,  1961). 

4.1.2  The  Instabilities  of  Wave  Breaking 

Great  progresses  have  been  made  in  recent  years  in  the  detailed 
analyses  of  periodic  deep  water  waves  of  their  amplitudes  up  to  the 
limiting  form,  their  responses  to  applied  surface  pressure  pulses,  their 
stability  properties  to  large  scale  perturbations,  the  structure  and  the 
stability  properties  of  their  small  scale  normal  modes,  and  their 
evolution  towards  wave  breaking.  The  whole  series  of  developments  on 
these  topics  are  based  upon  the  ability  to  calculate  the  free  surface 
deformation  (Longuet-Higgins  and  Cokelet,  1976)  of  an  irrotational 
periodic  motion. 

4. 1.2.1  Large  Scale  Subharmonics 

Longuet-Higgins  (1978b)  considered  a  set  of  M  number  of  waves  in  a 
spatial  interval  containing  NM  cycles  of  the  perturbation  waves.  He  was 
able  to  show  that  the  coalescence  of  the  frequencies  and  the  onset  of 
Benjamin-Feir  instability  are  delayed  to  larger  wave  slope  as  the  relative 
change  in  wave  number  increases.  One  interesting  result  is  that  once 
the  growth  rates  of  these  subharmonic  (large  scale)  waves  reach  the 
unstable  states,  they  will  restabilize  themselves  at  larger  value  of 
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wave  slopes.  Mathematically,  the  value  of  wave  slope  where  wave  breaking 
may  occur  has  been  suggested.  Consequently,  Longuet-Higgins  and  Cokelet 
(1978)  trace  the  initially  small  unstable  disturbances  upon  the  finite 
amplitude  wave  all  the  way,  ultimately,  to  breaking.  Their  results 
indicate  that  it  is  necessary  to  consider  that  the  short-time  behavior 
of  waves  and  the  surface  configurations  of  the  wave  crests  before 
breaking  are  similar  under  a  wide  range  of  physical  conditions.  The 
local  theory  of  wave  breaking  has  not  been  achieved,  however. 

4. 1.2. 2  Small  Scale  Superharmonlcs 

Longuet-Higgins  (1978a)  also  derived  the  normal  mode  perturbations 
of  a  nonlinear  gravity  wave  for  smaller  scales  than  that  of  the  basic 
wave.  As  the  steepness  of  the  basic  wave  increases  the  frequency  of 
each  perturbation  shifts  to  the  lower  value. 

4.1.3  Long  Wave-Short  Wave  Interactions 

With  the  conservative  Equations  (4.1)  and  (4.2),  encouraged  by 
Longuet-Higgins'  (1978a)  results  on  superharmonics,  the  interactions 
between  relative  long  and  short  waves  can  be  considered.  Of  course,  the 
accuracy  of  analytical  procedures  depends  asymptotically  on  the  ratlu  of 
long  to  short  wavelengths.  This  ratio  is  equivalent  to  the  mode  number 
N  defined  in  Section  4. 1.2.1. 

Considering  a  train  of  short  waves  interacting  with  the  underlining 
long  wave  in  a  Lagrangian  reference  frame  moving  with  the  long  wave,  the 
local  energy  density  of  the  snort  wave  train  is  simply 


E  =  \  P§a2 

(4.5) 

-  1  n^22 
-  -  pa  a 

(4.6) 

2|k| 

where,  neglecting  capillary  waves, 

a  =  (gj£|)1/2 


(4.7) 


I+1 2 

9  =  g  cos0  + 


(4.8) 


p  is  the  density,  a  is  the  wave  amplitude,  a  is  the  intrinsic  wave  frequency, 
j£|  is  the  wave  number,  0  is  the  angle  of  wave  slope  with  the  horizontal 
plane,  R  is  the  surface  curvature,  and  |u|  Is  the  tangential  speed  at  the 
surface  In  the  reference  frame.  From  Bernoulli's  Equation, 

|u(s)|2  =  |?|2  -  ?gc  (4.9) 

where  s  is  a  measure  of  distance  along  the  surface  of  long  wave,  |\*|  is  the 
phase  speed  of  the  long  w?.ve,  and  ?  is  the  free  surface  elevation. 

For  both  waves  travelling  in  the  same  direction,  Equation  (4.1)  gives 

a  +  |l<|  |u(s)|  =  am  +  |itj  |?|  =  aQ  (4.10) 

where  0Q  is  a  constant,  and  am  and  |l<m|  are  the  intrinsic  wave  frequency  and 
wave  number  at  the  mean  water  level.  From  Equation  (4.7),  |£|  =  o2/g, 
Equation  (4.10)  becomes 

( |u(s ) |/g)  a2  +  a  -  aQ  =  0 


or 


5(s)  =  (g/2|u(s) | ){-  1  +  [l  +  4a0|u(s)|/g]1/2}  (4.11) 

wnere  aQ  is  a  constant,  and  g  and  |u(s)|  are  functions  of  distance  s  along  the 
long  wave  free  surface.  The  intrinsic  wave  number  |k|  can  be  found,  subse¬ 
quently,  from  Equation  (4.10) 

|it(s)|  =  - (g/2ju(s)  I ){-  1  +  [1  +  4o0  |u(s)|/0]1/2}  (4.12) 

|u(s)| 

Comparing  the  values  of  |lt|  given  by  Equation  (4.12)  along  s  against  those 
calculated  by  Longuet-HIggins  (1978a)  for  small  scale  superharmonics,  it 
has  been  found  that  Equations  (4.11)  and  (4.12)  are  indeed  valid  for  the 
ratio  of  long  wavelength  to  short  wavelength  as  small  as  3  to  1  (Phillips, 
1979). 
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From  Equation  (4.2),  the  conservation  of  wave  action  flux,  the  short  wave 
energy  density  can  be  specified  as 


( |u(s)  I  ♦  (|f|  ^ 

y  a(s) 


(4.13) 


m 


where  |c  I  and  E„  are  the  phase  speed  and  the  short  wave  energy  density  at  the 
1  m1  m 

mean  water  level. 

From  Equations  (4.11),  (4.12),  (4.13)  and  (4.6)  for  d(s)/am, 

|*t(s)|/|lcm| ,  E(s)/Em,  and  (a|lc(s)  |  )2/(a|lcm|  )2  it  Is  important  to  notice  that 
small  change  In  long  wave  slope  can  produce  large  change  in  the  mean  square 
slope  of  short  wave. 


4.1.4  Interactions  Between  Short  Waves,  with  Spectral  Band,  and 
Long  Wave 

The  variations  of  the  apparent  wave  frequency  with  the  phase  of  the 
underlining  long  wave  observed  In  Eulerian  reference  frames  are  the  result  of 
two  effects.  One  is  the  variation  in  intrinsic  wave  frequency  o,  discussed  in 
Section  4.1.3,  and  the  other  is  the  variation  due  to  Doppler  effect  associated 
with  the  convection  of  short  waves  by  the  orbital  velocity  of  the  long  under¬ 
lining  wave.  The  apparent  wave  frequency  a  can  be  obtained  as 

a 


s  -  -  H  ■ 3  +  ^  fo +  i“(s>icos9) 


where  x  is  the  phase  function  of  the  long  wave  and  the  rest  of  the  symbols  are 
defined  as  before.  It  should  be  noted,  however,  |£(s)|/cos0  is  the  magnitude 
of  the  horizontal  component  of  |£(s)|  and  (-  |?|  ▼  |u(s)|cos6)  is  the  magni¬ 
tude  of  the  horizontal  component  of  the  long  wave  orbital  velocity  at  the 
surface.  From  Equation  (4.10) 


5/  ~ 

"  3t  '  0 


liUIL  |t| 

cose  11 


(4.14) 


The  last  term  at  the  right  hand  side  Is  simply  a  Doppler  shift  modification  on 
Equation  (4.10)  associated  with  the  change  of  reference  frame  from  that  of 
Lagrangian  to  Eulerian.  The  Lagrangian  reference  frame  Is  moving  at  the 
constant  phase  speed  |£|  of  the  long  wave. 
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For  this  short  wave  train,  at  a  fixed  position  x, 

?(x,t)  ■  a(x  -  ?  t)eix  (4.15) 

at  a  fixed  frequency  m  the  Fourier  component 

F(x,w)  ■  c(x,t)e1utdt 

•  }  /J  a(x  -  tt)e1+dt  (4'16) 

where  T  Is  the  long  wave  period  and  <j>  s  x  +  wt.  It  Is  quite  obvious,  from 
Equation  (4.14),  that  m  should  lie  at  and  between  the  extremes  of  -3x/3t. 
Phillips  (1979),  employing  the  stationary  phase  method,  obtains  the  spectral 
density  4(oi)  for  w  lying  between  the  extremes  of  -3x/3t  and  also  obtains  the 
spectral  densities  for  in  lying  at  the  extremes  of  -3x/3t.  The  spectral 
transformation  of  a  short  wave  train  in  the  presence  of  long  wave,  in  general, 
is  distributed  over  a  range  of  frequencies  of  width  (cr|lc|)|£|  where  6|t| 
is  the  difference  In  short  wave-numbers  between  at  the  crest  and  at  the  trough 
of  the  long  wave  and  |fj  Is  the  long  wave  phase  speed.  The  maximum  spectral 
density  occurs  at  the  highest  wave  frequency,  and  is  derived  for  short  waves 
at  the  crest  of  the  long  wave. 

4.1.5  Interactions  Between  Microscale  Surface  Waves  and  Long  Wave 

At  the  crest  of  the  long  wave.  It  is  Indicated,  In  the  previous  section, 
that  the  slopes  of  short  waves  will  increase.  If  these  wave  slopes  (for 
wavelengths  between  5  to  30  cm)  become  sufficiently  large,  the  short  waves  may 
splash  down,  produce  microscale  breaking,  or  may  form  parasitic  capillary 
waves . 

The  splashing  phenomenon  of  small  gravity  waves  has  been  observed  near 
the  crests  of  the  long  waves  and  cannot  be  explained  as  yet.  Microscale  wave 
breaking  requires  the  presence  of  the  wind-driven  current  (Banner  and  Phillips, 
1974;  PhllliDS  and  Banner,  1974)  and  can  be  enhanced  by  longer  gravity  waves 
somewhat.  Microscale  wave  breaking  can  also  be  suppressed  by  longer  gravity 
waves  when  the  wind  stress  Is  sufficiently  large  enough  to  regenerate  the 
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surface  wind-driven  layer.  The  thin  thermcl  conductive  layer  at  the  water 
surface  can  be  disrupted  by  microscale  wave  breaking.  The  formation  of 
parasitic  capillary  waves  ahead  of  short  gravity  wave  crests  has  been  analyzed 
theoretically  by  Longuet-HIgglns  (1963).  The  phenomenon  Is  conmunly  observed 
but  first  measured  by  Chang  et  al.  (1978)  In  a  laboratory.  There  are  qualita¬ 
tive  agreements  on  the  distribution  of  capillary  wavelengths  and  attenuation 
rate  with  respect  to  distance  along  the  long  wave  profile  just  ahead  of  the 
crest.  The  amplitudes  of  the  capillary  waves  were  measured  to  be  several 
orders  of  magnitude  larger  than  those  predicted  theoretically  by  Longuet- 
HIgglns. 

In  field  conditions  these  three  processes  may  all  occur  at  different 
proportions  with  various  wind  speeds.  Nevertheless,  their  spatial  densities 
are  the  Indicators  of  local  wave  energy  dissipation  at  capillary  scale  and, 
also,  short  gravity  waves  in  the  wavelengths  of  5  to  30  cm. 

4.2  Application  of  Long  Wave- Short  Wave  Interactions  to  Internal  Wave 

Modulations  of  Surface  Wave  Field 

In  Section  4.1.1  brief  mention  was  made  of  the  phenomenon  when  a  current 
distribution  moves  through  the  ambient  fluid  such  as  when  an  Internal  wave 
pulse  or  packet  moves  across  the  ocean  surface.  If  |£|  Is  the  speed  of  the 
internal  wave  pulse  (or  packet)  pattern,  the  surface  current  field  associated 
will  have  the  distribution  of  0  3  ft(x  -  |*|t)  along  the  x  direction.  The 
modulation  in  surface  wave  energy  density  E  Is  given  by  Equation  (4.4)  whore 
| c | ,  evaluated  from  Equation  (4.1),  is  the  phase  speed  of  the  surface  wave  and 
c,  ?,  and  ft  are  all  colllnear  in  the  direction  of  x  coordinate,  say.  For 
Internal  waves,  |0|  «  |c|  or  |tj.  As  the  group  speed  ^|c|  approaches  the 
propagation  speed  of  Internal  wave  pattern  |?|,  the  surface  wave  modulation 
will  increase  its  intensity.  This  physical  condition  was  analyzed  by  Phillips 
(1974)  in  which  the  situation  where  c,  t,  and  0  are  not  colllnear  has  also 
been  described.  Lewis  et  al.  (1974)  dar/vjstrated,  in  their  theoretical  and 
experimental  study,  the  existence  of  the  modulation  at  |C|  *  j|c|,  the  so- 
called  "Resonant"  or  "Blockage"  condition „  and  obtained  good  agreement  quanti¬ 
tatively  between  their  theory  and  exneriments. 
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Tha  velocity  ft  and  strain  rates  Induced  at  the  ocean  surface  by  the 
Internal  wave  are  small.  It  can  Induce^ surface  wave  modulation  significantly 
only  when  the  group  velocity  of  the  surface  wave  Is  fairly  close  or  equal  to 
the  propagation  speed  (or  phase  speed)  of  the  Internal  wave»  according  to  the 
resonant  (or  blockage)  condition.  The  fastest  phase  speed  of  the  Internal 
wave  (lowest  mode)  Is  of  the  order  of  50  cm/sec.  This  condition  certainly 
places  a  restriction  on  the  wavelengths  of  surface  waves  which  would  have 
significant  modulations  at  60  cm  or  less.  Generally,  the  ocean  surface  wave 
spectrum  contains  much  longer  wave  components  than  those  In  the  wavelengths  of 
60  cm  or  less.  Since  these  longer  waves  ar^  ore  energetic,  the  short  waves 
are  In  fact  strained  and  distorted  by  thes-  longer  waves  In  a  manner, 
described  In  Section  4.1.4,  far  more  rapidly  and  extensively  than  they  are  by 
the  Internal  wave  Induced  velocity  and  strain  rates.  Nevertheless,  the 
surface  wave  modulations  by  the  internal  wave  can  still  be  shown  observable. 

Let  the  surface  displacement  of  this  long  surface  wave  be  c  with  the 
phase  function  x  such  that 

;(x,t)  a  a  cosy, .  X  *  |£|x  -  crt  (4.17) 

which  Is  travelling  also  along  the  same  direction  x  as  the  internal  wave.  The 
short  gravity  wave  of  wave-number  lc'(s)  and  energy  density  E'(s)  Is  riding  on 
this  long  wave  in  a  celllnear  direction.  Both  £'( s)  and  E'(s)  of  tills  short 
gravity  wave  vary,  of  course,  with  respect  to  the  phase  function  x  of  the  long 
wave,  described  In  Section  4.1.3,  and  position  x  within  the  internal  wave 
current  pattern.  In  the  lagrangian  reference  frame  moving  with  the  internal 
wave  phase  speed  |ft|,  the  surface  current  induced  by  the  Internal  wave  is 
|ft|  -  |ft|  and  the  phase  function  x  is  changed  to  x  ■  |£|[x  -  (|c|  +  |ft|  -  |ft|)t] 
where  c  is  the  phase  velocity  of  the  long  wave.  According  to  Phillips  (1979), 
the  key  to  the  existence  of  an  analytical  solution  to  this  long-  short-internal 
wave  Interaction  problem  is  that  the  Intrinsic  wave  frequency  a"(s )  for  the 
short  gravity  waves  Is  constant  to  the  order  of  a | it|  of  the  long  wave  within 
Its  phase  function  cycle.  The  presence  of  long  wave,  therefore,  produces  no 
modulation  In  o'  at  this  order  for  the  short  gravity  waves.  The  significant 
modulation  in  o'  at  the  order  of  a|lc|  will  only  arise  from  the  short  wave- 
current  Interaction  or  o'  ■  o'(x).  Therefore,  the  surface  wave  modulations  by 
Internal  wave  Is  theoretically  observable. 
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5.  SURFACE  EFFECTS  DUE  TO  BATHYMETRY,  SUBSURFACE  MESOSCALE  FEATURES , 

INTERNAL  WAVES,  SUBSURFACE  EKMAN  LAYER  AND  BREAKING  WAVES,  AND 

THEIR  RELATIONSHIP  TO  THE  SUBSURFACE  PROCESSES 

In  the  previous  three  chapters,  the  responses  of  the  ocean  surface  waves 
to  a  moving  current  pattern  or  stationary  current  variations  (or  gradients) 
generated  by  subsurface  dynamic  processes  have  been  discussed  briefly.  There 
are  two  areas,  pertinent  to  this  survey,  required  to  be  addressed. 

The  first  area,  of  course.  Is  what  kind  of  accuracy  can  be  expected  on 
the  surface  current  variations  derived  from  the  observed  surface  wave  effects. 
Intuitively,  this  accuracy  Is  of  the  order  of  a | tc |  because  the  surface  current 
variations  would  have  to  be  derived  through  strong  interaction  phenomena. 

The  second  area  of  concern  Is  what  can  be  achieved,  by  knowing  the 
surface  wave  and  the  surface  current  variation  fields,  in  resolving  the  sub¬ 
surface  features  such  as  bathymetry,  mesoscale  subsurface  eddies  and  fronts, 
internal  waves,  and  subsurface  shear  currents,  for  example.  The  validities  of 
the  inverse  algorithms  which  would  have  to  be  used  to  Infer  the  subsurface 
features  depend  on  how  complete  the  mathematical  field  equations  would  be  for 
modelling  the  dynamic  physical  phenomena.  The  complete  set  of  mathematical 
field  equations  which  are  briefly  summarized  by  Chen  and  Noble  (1980)  govern 
the  three-dimensional  current  and  density  fields  in  the  ocean,  and  the 
kinematic  and  the  dynamic  boundary  conditions  at  the  oceanic  boundary.  These 
boundary  conditions  describe  the  states  of  heat,  mass,  and  momentum  exchanges 
as  functions  of  position  and  time.  At  the  ocean  surface  these  surface  boundary 
conditions  include  strong  surface  wave  interaction  phenomena  to  the  order  of 
a |lc |  and  weak  surface  wave  interaction  phenomena  to  the  order  of  (a|lc|)2. 

Heat  exchange  surface  effects  (Phillips,  1974)  are  weak  interaction  phenomena 
while  mass  exchange  surface  effects  can  be  strong  (e.g.,  estuarine  or  coastal 
run-off)  or  weak  (e.g.,  evaporation)  interaction  phenomena. 

The  steady  uniform  current  field  creates  no  surface  wave  effects  but  has 
dynamic  topographic  signatures  (Hill,  1962)  through  the  geostrophic  assump¬ 
tions  (baroclinlc  and  barotroplc).  Various  types  of  dynamic  topography  have 
been  observed  over  currents  (Chen  and  Noble,  1980)  and  subsurface  mesoscale 
eddies  and  fronts  (Chen  and  Noble,  1980;  Dugan  and  Chen,  1980).  The  treatise 
on  dynamic  topography  for  the  Inference  of  geostrophic  features  such  as 
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currents  and  eddies  is  discussed  in  detail  by  Chen  and  Noble  (1980).  In  a 
special  situation  where  barotroplc  Instability  Is  assumed  to  be  excited  by 
tidal  oscillations,  Mollo-Chrlstensen  et  al.  (1981)  successfully  give 
estimates  of  current  speed  In  the  Gulf  Stream. 

5.1  Bathymetry 

As  ocean  surface  waves  propagate  from  deep  to  shallow  water  with  variable 
depth,  significant  changes  occur  In  wave  spectral  characteristics.  These 
changes  are  the  combined  effects  due  to  physical  dynamic  processes  which  may 
or  may  not  be  governed  by  the  conservation  of  energy.  The  physical  dynamic 
processes  whose  total  wave  energy  densities  are  conserved  Include  those 
processes  such  as  wave  refraction,  wave  shoaling,  and  wave-wave  Interaction. 
The  physical  dynamic  processes  whose  total  wave  energy  densities  are  net 
conserved  can  be  those  Involved  in  energy  dissipation  and  generation. 

5.1.1  Wave  Refraction  and  Shoaling 

Based  on  the  conservation  of  energy  flux,  Longuet-Higgins  (1957)  consi¬ 
dered  the  transformation  of  continuous  wave  spectrum  by  analysis  of  the 
refraction  process  for  just  two  components.  As  an  extension,  Karlesson* 

(1969)  generalized  the  transformation  for  the  continuous  directional  wave 
spectrum  by  nalysls  of  its  distribution  under  steady  state  conditions  with 
finite  difference  techniques.  Using  a  geometrical -optical  approximation, 
Krasltskiy  (1974)  derived  the  explicit  analytical  solution  for  the  spectral 
transformation  over  two-dimensional  parallel  bottom  contours. 

The  laboratory  and  field  experiments  which  were  designated  to  observe  the 
phenomena  in  wave  spjetra!  transformation  include  Wu  et  al.  (1977)  and  Shemdln 
et  al.  (i9oC)  among  others.  Bear  In  nv»r.d  that  the  data  observed  contains 
information  on  phenomena  beyond  wave  refraction  and  shoaling.  Wu  et  al.  (1977) 
made  come  laboratory  observations  on  wave  energy  dissipation,  wave  energy 
generation,  and  also  on  nonlinear  waye-wave  Interaction  phenomena  in  the 

*Collins  (1972)  further  extended  the  transformation  to  include  energy 
dissipation  and  generation  and  applied  his  results  to  an  Irregular  bottom 
bathymetry.  By  tracing  wave  energy  along  wave  rays,  his  numerical  scheme, 
however,  made  the  computational  procedures  quite  Impractical  for  the 
determination  of  shallow  water  wave  spectrum  at  predetermined  and  designated 
locations. 
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context  of  wave  spectral  transformation.  Recently,  Shemdln  et  al.  (1980) 
summarized  some  of  the  linear  and  nonlinear  physical  dynamic  processes  which 
a**e  relevant  to  the  wave  spectral  transformation  In  finite-depth  water  and 
applied  the  results  against  data  obtained  from  the  JONSWAP  (Hasselmann  et 
al.,  1973)  site  In  the  North  Sea.  They  concluded  that  the  rate  of  wave  energy 
dissipation  Is  dominated  by  yet  unidentified  physical  dynamic  processes  which 
are  determined  by  the  bottom  conditions. 


5.1.2  Wave  Energy  Dissipation* 

There  are  five  different  physical  dynamic  processes  associated  with  wave 
energy  dissipation.  They  are:  Percolation,  Wave-Induced  Bottom  Motion,  Wave 
Attentuatlon,  Bottom  Friction,  and  White-caps  and  Wave  Breaking.  All  but  the 
last  two  are  linear  processes  (or  mechanisms). 


5. 1.2.1  Percolation 


The  wave  pressure  field  induces  a  flow  In  the  sandy  bottom  (or  bed) 
below.  The  magnitude  of  this  induced  flow  is  a  dependent  of  the  permeability 
of  the  sandy  bed.  The  wave  energy  dissipation  can  be  computed  from  the  cross¬ 
correlation  between  the  wave  Induced  pressure  above  the  bed  and  the  vertical 
velocity  just  at  the  top  of  the  bed.  Putnam  (1949)  made  the  first  investi¬ 
gation  for  an  isotropic  sand  layer.  McClain  et  al.  (1977)  made  corrections  to 
the  model  regarding  the  continuity  at  the  bottom  interface.  Hsiao  and  Shemdln 
(1978)  improved  Putnam's  results  to  accommodate  different  permeabilities  In 
the  horizontal  and  the  vertical  directions.  The  rate  of  energy  dissipation 
for  water  depth  h  Is  given  by 


3E  _ 


3t  E | k | 


tanh/(<H/<y)  |  k  J  d 


coshz|l?|h 


(5.1) 


where  icH  and  <v  are  the  horizontal  and  vertical  permeability  coefficients 
(Darcy's  numbers),  respectively,  and  d  is  the  thickness  of  the  sand  layer. 


*Wave  Energy  Generation  Is  beyond  the  scope  of  this  survey. 
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5. 1.2. 2  Wave- Induced  Cotton  Motion 

The  bottom  responds  to  the  wave-induced  pressure  field  with  viscoelastic 
characteristics  when  the  material  of  bottom  sediment  Is  mud  or  decomposed 
organic  matters.  Some  field  and  laboratory  observations  have  been  made  on 
this  phenomenon  by  e.g.,  Tubman  and  Suhayda  (1976),  Rosenthal  (1978)  among 
others.  A  solution  for  a  two-layer  coupled-flow  model  was  given  by  Hsiao  and 
Shemdln  (1980)  for  a  finite  depth  mud  layer. 

5. 1.2. 3  Wave  Attenuation 

The  direct  attenuation  by  molecular  viscosity  strongly  Influences  the 
state  of  capillary  wave.  If  a  slick  covers  the  surface,  the  amplitude  of 
shorter  gravity  waves  may  also  be  attenuated  as  well.  On  the  longer  gravity 
waves,  the  direct  effect  of  molecular  viscosity  is  negligible.  Phillips  (1959) 
Investigated  this  physical  mechanism  extensively. 

5. 1.2. 4  Bottom  Friction 

Energy  dissipation  due  to  bottom  friction  primarily  Is  the  result  of  two 
physical  mechanisms.  One  Involves  the  dissipative  work  done  against  turbulent 
shear  stresses,  the  so-called  turbulent  bottom  friction,  induced  at  the  bottom 
bed  by  water  particle  motions.  The  other  Is  responsible  for  the  dissipative 
work  done  against  the  viscous  forces  induced  at  the  permeable  bottom  where 
percolation  occurs. 

In  general,  turbulent  bottom  friction  Is  dominant  over  other  dissipative 
processes  when  the  sediment  Is  composed  of  sand  with  mean  diameter  In  the 
range  of  0.1  -  0.4  mm  where  low  permeability  prohibits  percolation  and  granular 
friction  prevents  viscous  flow  behavior.  Putnam  and  Johnson  (1949)  Investigated 
this  problem  Initially.  They  used  the  quadratic  friction  law  to  derive  the 
rate  of  energy  dissipation  for  sinusoidal  surface  waves.  Assuming  Gauss Un¬ 
distributed  surface  wave  field,  Hasselmann  and  Collins  (1968)  derived  the  rate 
of  energy  dissipation  for  a  random  sea.  With  the  concern  placed  on  bottom 
friction  coefficient,  Bretschnelder  and  Reid  (1954)  evaluated  the  coefficients 
for  the  sandy  bottoms.  Furthermore,  relying  on  labontory  and  semi-theoretical 
studies,  Johnson  (1965)  found  that  bottom  coefficients  are  functions  of  Reynolds' 
number,  and  also  of  relative  roughness.  Hsiao  and  Shemdln  (1978)  have  Investi¬ 
gated  this  problem  In  detail. 


19 


5. 1.2. 5  White-Caps  and  Wave  Breaking 

It  Is  a  well-observed  fact  that  under  sustained  wind  action,  waves 
grow  to  some  maximum  amplitudes  before  the  onset  of  breaking  or  white- 
caps  would  limit  their  growth.  Hasselmann  (1974)  treated  the  white-cap 
Interactions  In  terms  of  an  equivalent  ensemble  of  random  pressure 
pulses.  By  assuming  Irrotatlonal  motion  Longuet-Hlgglns  and  Cokelet 
(1976;  1978)  calculated  the  deformation  of  the  free  wave  surface  up  to 
breaking  for  horizontally  periodic  surface  waves.  Subsequently,  Longuet- 
Hlgglns  (1978a;  1978b)  calculated,  for  the  periodic  deep  water  waves, 
their  amplitudes  up  to  limiting  forms,  their  responses  to  applied  surface 
pressure  pulses,  their  stability  properties  to  subharmonic  (large  scale) 
perturbations,  the  structure  and  stability  of  their  superharmonic  normal 
modes  and  their  evolution  toward  breaking.  The  results  provide  great 
progress  In  this  field.  The  subjects  were  covered  more  extensively  In 
Section  4.1.2  In  the  chapter  on  strong  Interactions.  Nevertheless,  more 
research  work  Is  still  required,  especially  for  surface  waves  over 
variable  depth. 

5.1,3  Wave-Wave  Energy  Transfer  In  Finite-Depth  Water 

The  nonlinear  wave  energy  transfer  due  to  weak  nonlinear  wave-wave 
Interactions  is  derived  from  Equations  (3.1),  (3.2)  and  (3.4)  and  Is 
given  by  a  Boltzman  Integral  of  the  form. 


(5.2) 

•  <$(l<i  +  ^2  "  ^3  "  Mai  +  *  °3  "  ) <l^i  die gdTc 3 


3F(k4,h) 


at 


hW*  +  °3F1F2F4  '  01F2F3F4  '  a2FlF3F4] 


where  the  rate  of  change  of  the  spectrum  at  wave-number  $4  due  to  inter¬ 
actions  with  other  components  lt2,  and  lt3  was  described  in  chapter  3. 

F^  Is  the  spectral  density  at  1T^,  or  F^  *  F(^).  The  depth  dependence 
in  Equation  (5.2)  enters  via  the  coefficient  D  and  the  relationship, 
given  by  Equation  (3.2)  and  (3.3),  particularly  In  the  frequency  delta 
function.  Exact  computations  of  this  Integral  were  performed  by  Hasselmann 
and  Hasselmann  (1980)  In  finite-depth  water  for  directional  wave  spectra. 
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In  the  presence  of  currents  the  radiative  transfer  equation  for  the 
wave  spectrum  F{x,  It,  t)  has  the  form. 


a 

■5 T 


<?> 


+  c. 


(5.3) 


where  cg  is  the  group  velocity,  a  is  the  Intrinsic  wave  frequency  given 
In  Equation  (4.1)  with  respect  to  the  reference  frame  moving  with  the 
current,  and  S  Is  the  net  source  function.  The  source  function  S  is  the 
sum  of  several  components  such  as 


S 


$.  *  S 
f  w 


(5.4) 


where  Sp,  Sfa,  are  associated  with  energy  dissipation  by  percolation, 

wave-induced  bottom  motion,  and  bottom  friction,  respectively.  Sw  Is 

associated  with  the  energy  transfer  by  weak  nonlinear  wave-wave 

Interaction.  S.  Is  associated  with  •.  e  attenuation,  and  S„  is  associated 
a  g 

with  energy  generation  by  wind. 

The  terms  F/a  is  changed  to  F  in  the  absence  of  current  systems. 

In  the  case  of  steady  wind-driven  and  thermohaline-driven  barotroplc 
velocity  or  current  field  In  a  variable-depth  ocean  under  the  mild 
assumption  that  the  bathymetry  does  not  penetrate  up  to  the  depth  of 
appreciable  horizontal  density  gradients,  the  bottom  torque  (or  vortlclty) 
can  be  related  directly  to  the  stretching  term  In  the  barotroplc  velocity 
balance  (Rattray  and  Dworskl,  1978).  Thei  Tore,  even  this  barotroplc 
current  field  can  be  modified  by  the  bathymetry  at  the  ocean  surface. 

The  plausibility  of  zero  current  field  description  over  a  coastal  region 
with  finite-depth  is  questionable. 


5.2  Subsurface  Mesoscale  Features 

Although  oceanic  features  such  as  currents,  mesoscale  eddies,  and 
mesoscale  fronts  are  shown  to  exhibit  measurable  temperature  gradients 
at  both  the  surface  and  subsurface  levels  (for  example,  Stommel,  1966; 
Stommel  and  Yoshida,  1972;  Bruce,  1979;  Halllwell  and  Mooers,  1979; 
Richardson  et  al.,  1979;  Vukovlch  and  Crlssman,  1978  ),  in  large  mid¬ 
ocean  areas,  poor  correlation  coefficients  of  the  order  of  0.2-0. 3 
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between  surface  temperature  patterns  and  the  thermocllne  structure  were 
found  by  Dugan  (1980)  and  Dugan  and  Chen  (1980).  A  majority  of  these 
oceanic  features  do  possess  measurable  dynamic  heights  and  slopes  (for 
example,  Stommel,  1966;  Stommel  and  Yoshida,  1972;  Vukovlch  and  Crissman, 
1978;  Dugan,  1980;  Dugan  and  Chen,  1980).  These  dynamic  heights  and 
slopes  are  the  result  of  different  density  and  current  profiles  in  the 
horizontal  and  vertical  directions  through  the  currents,  the  eddies,  and 
to  some  extent  the  fronts.  These  dynamic  heights  and  slopes  are  the 
deviations  from  the  local  geold  heights  and  slopes  which,  theoretically, 
represent  a  motionless  ocean.  The  dynamic  slopes  as  well  as  the  variations 
of  dynamic  heights  can  induce  horizontal  pressure  gradients.  Under  the 
steady  state  situation,  by  neglecting  the  effects  due  to  friction  and 
removing  the  effects  due  to  wind  set-up  and  atmosphere  pressure,  the 
geostrophic  current  velocity  is  related  orthogonally  to  the  horizontal 
pressure  gradient  by  the  dynamic  method.  In  the  cac«  of  eddies,  the 
horizontal  pressure  gradient  is  further  balanced  by  the  centrifugal 
force.  Therefore,  under  the  geostrophic  approximation,  the  current 
velocity  field  can  be  evaluated  if  the  dynamic  heights  and  slopes  over 
the  current,  the  eddy,  or  the  fronts  are  known. 

Theoretically,  a  complete  sec  of  conditions  required  for  a  mathe¬ 
matically  and  physically  unique  and  bounded  current  velocity  field  are 
well  defined.  These  conditions,  derived  from  the  Euler's  equation  of 
motion  and  the  equation  for  continuity  of  fluids,  are  the  known  bottom 
topography,  the  known  vertical  and  horizontal  density  profiles,  and  the 
known  boundary  kinematic  and  dynamic  conditions.  In  other  words,  if  the 
current  velocity  field  is  known,  the  unknown  three-dimensional  density 
field,  which  Includes  the  subsurface  mesoscale  eddies  and  fronts,  can  be 
uniquely  determined. 

Assuming  a  linear  constitutive  law  for  the  surface  strain  rate 
tensor,  the  basic  hydrodynamic  equations  describing  oceanic  circulation 
in  a  general  time  dependent  case  are  (Hill,  1962)  the  Euler's  equation 
of  motion, 
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(5.5) 


[0(x,t)*v]0(x,t)  +  20x0(x,t)  +  w(x,t)x(w(x,t)xx] 

* - -r —  0p(x,t)  +  0(x,t)  +  vV20(x,t)  * 

p(x,t) 

and  the  equation  of  continuity, 

+  |>(x.t)ff(x.t)]  =  0  (5.6) 

with  the  associated  kinematic  and  dynamic  boundary  conditions.  In  Equations 

(5.5)  and  (5.6),  0  is  the  current  velocity,  $  is  the  angular  velocity  of 
the  earth,  u>  is  the  angular  velocity  of  the  eddy,  x  is  the  position 
vector,  p  is  the  density  of  sea  water,  p  is  the  pressure,  &  is  the  body 
force,  and  v  is  the  kinematic  viscosity  of  sea  water.  The  position 
vector  x  is  a  function  of  time.  The  third  term  on  the  left  hand  side, 

20x0,  of  Equation  (5.5)  is  the  Coriolis  force  induced  by  the  fluid 
motion  0  and  the  fourth  term,  uix(Jxx),  is  the  centrifugal  force  induced 
by  the  eddy  motion.  In  the  absence  of  eddy  motion  the  centrifugal  force 
term  always  vanishes.  For  steady  state  situation.  Equations  (5.5)  and 

(5.6)  become 

(0«^)0  +  20*0  +  (w«x)w  -  (u!*u>)x  =  -  1  0p  +  d  +  vV20  (5.7) 


and 


0*0p  +  p^*0  =  o 


(5.8) 


With  the  Boussinesq,  the  hydrostatic,  the  Rectangular  Cartesian  Coordinates', 
and  the  scale  approximations  and,  also,  the  eddy  term  assumption,  Equations 
15  7)  and  (5.8)  are  further  reduced  to  the  component  form  as 
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(5.10) 
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au2(x) 

+  3x2  " 


0, 


(5.12) 


respectively,  where  x  *  (x-j,  x2,  x3),  A  is  the  latitude  (positive  In  northern 
hemisphere)  and  g  is  the  gravitational  acceleration.  The  Cartesian  Coordinates 
x-|,  x2»  and  x3  are  In  the  local  east,  north,  and  vertical  downward  directions. 
However,  since  usually 


2  =  (0,0,u>3)  (5.13) 

Equations  (5.9)  to  (5.12)  can  be  simplified  as 


2sinA|&|U-(x)  +  w-A,  -  — ^ - 

C  p(x) 

.Mi). 

3x-j 

(5.14) 

2sinA|^|U, (x)  -  uj,2x0  = - — 

1  3  2  P(x) 
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3x2 

(5.15) 

^  *  p(*)« 

(5.16) 

3iyx)  3U2(x) 
3x-j  h  3x2 

a  0 

(5.17) 

Equation  (5.16)  is  the  hydrostatic  equation.  In  the  case  of  eddy  motion. 
Equations  (5.14)  to  (5.17)  become 

2s1nA|ft||ft(r)|  +  MrllL  »1  |£  (5.18) 
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(5.19) 


p  H  ■  9  (5.20) 

and 

.^(rlL  ,  o  (5.21) 

where  r,  0,  and  z  are  the  coordinates  for  the  cylindrical  coordinates  along 
the  radial,  angular,  and  vertical  directions. 

The  second  terms  at  the  1=ft  hand  sides  of  Equations  (5.14)  and  (5.15) 
are  high  order  terms,  for  current  systems  other  than  mesoscale  eddies,  these 
terms  can  be  dropped  as 

2s1nX|$|U?(x)  *  —V  ^T1  (5-22) 

2  P(x)  3xl 

2sinA|j$|U,  (x)  =  -  •  (5.23) 

1  P(x)  Sx2 


Equations  (5.22),  (5.23),  (5.16),  and  (5.17)  are  the  governing  equations 
for  the  geostrophic  current  systems.  From  Equation  (5.16)  or  (5.20),  let 
Equation  (5.16)  be 

P(x)  -  pa(x)  =  g  /x3  p(x)dx3  (5.24) 

where  p,(x)  is  the  reference  atmospheric  pressure  at  the  mean  sea  surface 

a 

S.  The  mean  sea  surface  s  Is  obtained  by  referring  the  dynamic  height  to 
an  arbitrary  datum  where  the  horizontal  pressure  gradient  and  current 
ve’ocity  are  zero  (Chen  and  Noble,  1980). 

Recent  results  of  numerical  multi-layer  model  studies  (Placsek,  1980) 
Indicate  that  the  subsurface  mesoscale  features  such  as  eddy,  a  predominant 
geostrophic  barocllnlc  feature,  should  exhibit  a  surface  effect  (or 
signature)  In  the  form  of  deviation  (ctynamlc  topography)  from  the  local 
geoid.  SEASAT-1  altimeter  has  been  able  to  observe  this  deviation  induced 
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by  the  edcly  1r.  the  Gulf  of  Mexico.  Through  the  surface  effects  due  to 
wave-current  interactions  (strong  Interactions)  the  surface  current  boundary 
conditions  of  these  subsurface  mesoscale  features  may  further  enhance  the 
possibility  of  resolving  their  three-dimensional  density  structures  to 
the  first  Internal  wave  mode  if  the  local  stability  (Brunt-VSIsSla} 
frequency  history  is  known. 


5.3  Internal  Wave  with  Subsurface  Shear  Current,  Layered  Density 
Anomaly,  Bathymetry,  and  Heat  Storage 

To  the  Bousslnesq  approximation,  and  neglecting  the  earth's  rotation, 
the  equation  describes  the  motion  of  a  stratified  fluid  having  the  form 
(Phillips,  1969)  of 

,2 


%  (v2u3)  +  N2(x3)72U3  *  Q(x,t) 
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where  the  nonlinear  terms, 
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(5.26) 


the  Brunt-Vais31a  frequency. 


the  continuity  equation, 
3IL 

a  n 

and 
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where 
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B  »  •  9  (5.29) 

o 

Is  the  reference  water  density  at  x3  «  0  or  at  the  ocean  surface. 
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For  Infinitesimal  disturbances  In  the  absence  of  a  mean  shear,  the  term 
Q(x,t)  in  Equation  (5.25)  Is  negligible  and  the  governing  equation  Is  linear; 
but  the  thermocline  Is  often  a  region  with  appreciable  mean  vertical  shear. 
Vertical  shear  has  the  effect  of  reducing  the  frequency  and  vertical  scale  of 
small  scale  internal  wave  motions  (Phillips,  1969).  In  the  absence  of  energy 
dissipation  phenomena,  energy  Is  eventually  lost  from  the  Internal  waves  and 
Is  being  transferred  to  the  mean  shear.  The  vertical  velocity  fluctuations 
are  thus  suppressed  and  the  motion  eventually  becomes  totally  horizontal. 

Bell  (1979)  used  real  background  oceanographic  data  and  examined,  by  numerically 
simulated  Internal  wave  patterns,  how  the  enerqy  envelope  (that  region  of 
space  within  which  appreciable  density  perturbations  occur)  and  the  phase 
Information  (the  structure  of  the  density  perturbation  within  the  energy 
envelope)  observed  In  the  absense  of  shear  are  distorted  by  the  Introduction 
of  background  «hear.  The  relevant  parameter  Bell  used  for  the  determination 
of  relative  importance  of  the  background  shear  is  the  bulk  Richardson  Number 
R(xg)  as 


R(x3) 


N2(x3) 

(3|tJh|/3x3)2 


(5.30) 


where  *  (U-|,U2).  Because  of  the  Influence  of  the  large  scale  averaged 
properties  of  the  shear  profile  on  the  wave  propagation  process.  Bell  found 
that  the  net  drift  velocity  of  the  energy  envelope  at  any  depth  level  does  no: 
accurately  reflect  the  background  current  speed  at  that  level.  Shear  Induced 
restrictions  of  the  vertical  extent  of  the  energy  envelope  are  shown  to  be 
consistent  with  the  wave  propagation  process  relating  to  the  critical  layer 
nhenomenon.  Significant  coherent  losses  relative  to  unsheared  patterns  over 
several  Brunt-Vaisala  periods  were  observed  on  the  phase  Information.  Further 
more,  the  energy  envelope  hardly  reached  the  ocean  surface. 

Based  upon  Bell's  numerical  simulation,  there  Is  no  surface  effect  on 
Internal  wave  processes  due  to  the  presence  of  the  subsurface  shear  current 
structures  or  patterns  except  that  the  amplitude  and  the  frequency  of  the 
Internal  wave  may  be  reduced  somewhat  through  the  high  order  mechanism.  It 
would  be  rather  difficult  to  observe  these  changes  In  amplitude  and  frequency 
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of  the  Internal  wave,  Induced  by  the  presence  of  subsurface  shear  current, 
through  the  mechanism  Identified  in  Section  4.2  above. 

An  Internal  wave  can  be  reflected  when  It  incidents  upon  a  layered 
density  anomaly  embedded  In  a  stably  stratified  fluid  with  constant  Brunt- 
VSIsalS  frequency.  Mled  and  Ougan  (1974)  found  that  the  reflection  coefficient 
Is  a  complicated  function  of  the  horizontal  and  vertical  wave  numbers  and  the 
parameters  related  to  the  Brunt-VSIsHlS  profile.  All  wave  energy  Is  trans¬ 
mitted  through  th*  anomalous  layer  for  waves  propagating  almost  horizontally 
and  a  significant  amount  of  wave  energy  Is  reflected  for  vertically  propa¬ 
gating  ones.  Through  WKB  description,  there  can  exist  for  some  stratification 
conditions  certain  other  directions  of  propagation  for  which  the  conditions 
for  energy  transfer  are  more  favorable  than  they  are  for  nearly  contiguous 
Incident  angles  of  Inclination  with  respect  to  the  horizontal.  The  existence 
of  these  reflected  Internal  waves  can  be  an  Indicator  of  the  presence  of  a 
layered  density  anomaly. 

Internal  waves  do,  of  course,  reflect  from  a  sloping  bottom  but  the 
significance  of  this  mechanism  to  ocean  surface  phenomena  is  much  less  Impor¬ 
tant  than  those  described  In  Section  5.1.  Internal  waves  can  also  be 
reflected  (Mollo-Chrlstensen,  1977)  by  a  moving  surface  layer  which  can  be 
Induced  by  spatially  varying  bathymetry,  for  example. 

In  the  Instances  where  the  Internal  waves  are  generated  by  tidal  flow 
over  bottom  topography  one  can  find  the  propagation  speed  of  these  Internal 
waves  from  their  wave  packet  spacing  (Hollo-Christensen  and  Mascarenhas, 

1979).  The  propagation  speed  depends  upon  the  density  anomaly  and  depth  of 
the  upper  mixed  layer.  Attributing  the  density  anomaly  to  temperature  rather 
than  to  salinity,  Mollo-Chrlstensen  and  Mascarenhas  can  calculate  the  heat 
storage  In  the  upper  oceanic  layer.  In  this  way,  heat  storage  was  calculated 
from  LANDSAT  satellite  data  over  the  New  England  continental  shelf  and  was 
compared  with  available  In-sltu  observations.  They  found  that  the  method  may 
have  merit  and  is  deservlno  of  further  refinement.  Furthermore,  If  one  could 
detect  surface  signs  of  Internal  waves  on  the  main  thermocllne  and  then  find 
that  such  waves  are  emitted  over  bottom  topography,  It  should  then  be  possible 
to  follow  the  integrated  density  anomaly  over  large  regions  of  the  ocean 
(Mollo-Chrlstensen,  1981). 
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5.4  Horizontal  Distortion  In  the  Surface  Ekman  Layer 

As  pointed  out  by  Mollo-Chrlstensen  (1981),  the  surface  wind-driven  Ekman 
layer  will  show  signs  of  variations  of  Internal  vortlclty  through  horizontal 
divergence,  and  signs  of  bottom  Ekman  layer  and  free  Internal  shear  layer 
variations.  Faller  and  Kay lor  (1966)  showed  that  Ekman  layers  develop  rolls; 
such  rolls  may  show  In  the  surface  Ekman  layer  as  Langmuir  circulations. 
Mollo-Chrlstensen  observed,  from  LANDSAT  Images  and  SEASAT-1  Synthetic  Aperture 
Radar  (SAR)  Imagery,  that  the  rolls  or  the  strlatlons  seem  to  be  parallel  to 
the  velocity  field  beneath  the  Ekman  layer,  rather  than  at  an  angle  with  the 
geostrophlc  flow  as  found  by  Faller  and  Kaylor  (1966).  These  strlatlons  may 
be  due  to  a  periodic  structure,  speculated  by  Mollo-Chrlstensen,  In  the  Ekman 
layers  at  the  bottom  or  the  Ekman  layers  at  shear  Interfaces  beneath  the 
surface. 

Mollo-Chrlstensen  does  not  have  an  analytical  solution  for  the  Ekman 
layer  created  by  stratified  flow  parallel  to  bottom  topography,  but  he  feels 
that  there  has  to  be  a  periodic  solution  because  of  the  effects  of  buoyancy 
that  limit  vertical  excursions  of  fluid  particles.  The  direction  of  Isobaths 
at  the  bottom  or  of  Isobaths  in  the  free  shear  layer  that  exist  below  the 
surface  current  can  be  deduced. 

5.5  Influence  of  Breaking  Waves  on  the-  Dynamics  of  the  Upper  Ocean 

Using  the  statistical  model  proposed  by  Longuet-Hlgglns  (1969),  Huang  et 
al.  (1981)  found  that,  based  upon  the  balance  of  energy  between  wind  and 
waves,  the  influence  of  breaking  waves  on  the  dynamics  of  the  upper  ocean 
depends  only  on  one  parameter,  the  significant  slope,  defined  and  used  by 
Huang  and  Long  (1980).  If  the  energy  released  from  the  breaking  waves  Is 
treated  as  the  sole  source  of  turbulence  energy,  various  models  can  be 
proposed  to  simulate  a  wide  variety  of  both  surface  and  subsurface  dynamic 
phenomena.  The  quantities  calculated  include  wave  attenuation  rate,  equi¬ 
valent  eddy  viscosity,  wave  Induced  drag  coefficient,  surface  drift,  white¬ 
capping  percentage,  and  mixing  efficiency  for  the  mixing  layer  models. 

The  treatment,  by  Huang  et  al.  (1981),  cf  using  breaking  waves  as  the 
turbulence  energy  source  alone  may  be  perceived  only  as  the  first  order  of 
approximation  of  the  quantities  calculated.  Therefore,  their  results  are  the 
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lower  bound  of  the  real  values  of  these  quantities.  The  quantitative  values 
of  all  their  results  are  found  to  agree  with  the  conventional  adapted  values 
very  well. 

The  Important  and  revolutionary  advantage  of  their  approach  Is  that  the 
Inputs  to  the  models  are  all  obtainable  from  remote  sensors.  Therefore,  the 
practical  implications  of  their  results  can  be  enormous  and  their  approach  Is 
deserving  of  further  theoretical,  experimental  and  field  Investigations. 

6.  GENERAL  SUMARY 

6.1  Weak  Interactions 

Two  types  of  phenomena,  based  upon  their  spectral  spreads,  can  be  catego¬ 
rically  classified  among  surface  wave  weak  interactions.  The  first  type, 
broad  band  phenomena.  Includes  wave-wave  energy  transfer  and  general  wave 
strain  Instability.  The  second  type,  narrow  band  phenomena,  includes  Benjamin- 
Felr  instability,  recurrence,  and  envelope  solitons.  All  these  phenomena 
share  the  physical  characteristics  of  the  following:  (a)  they  may  be  manifested 
only  by  the  long  dominant  waves  In  a  wave  field  and  yet  the  behavior  of  short 
waves  is  determined  by  the  strong  interactions  with  the  long  waves  rather  than 
by  these  processes;  (b)  they  are  evolutionary  phenomena  with  a  time  scale  of 
(ak)  times  the  wave  period,  where  a  Is  the  wave  amplitude  and  k  is  the  wave 
number  of  the  dominant  wave.  In  a  typical  oceanic  condition,  this  time  scale 
is  roughly  of  the  order  of  at  least  100  wave  periods  of  the  dominant  wave. 
Equivalently  speaking,  the  space  scale  is  of  the  order  of  at  least  100  wavelengths 
of  the  dominant  wave.  Therefore,  due  to  their  modifications  to  wave  field, 
which  are  of  the  same  order  of  (ak)"2,  weak  interactions  cannot  be  used  to 
account  for  any  local  disturbance  to  the  wave  field. 

6.2  Strong  Interactions 

The  strong  Interactions  are  those  phenomena  whose  time  scale  is  of  the 
order  of  the  wave  period  and  space  scale  Is  of  the  order  of  the  wavelength. 

The  detailed  wave  structures  (or  profiles)  are  the  prime  sources  of  Information. 
Therefore,  the  techniques  of  Fourier  Transform  representation  of  wave  field, 
mathematically,  is  not  valid  in  the  interpretation  for  strong  interactions. 


Strong  Interactions  Include  the  strong  Longuet-HIggins  Instability,  wave- 
current  Interactions,  long  wave-short  wave  Interactions,  and  the  processes  of 
parasitic  capillary  formation  and  microscale  breaking  Induced  by  surface  wind 
drift.  Due  to  the  rapid  surface  wave  responses  to  local  oceanic  disturbances 
accounted  by  the  strong  Interactions,  their  phenomena  can  offer  at  least 
qualitatively  as  the  Indicator  of,  for  example,  the  existence  of  Internal  wave 
patterns.  Additional  research  work,  absolutely  necessary  for  the  better 
quantification  of  these  phenomena,  Is  currently  hindered  by  the  lack  of  mathe¬ 
matical  tools.  Therefore,  quantitative  laboratory  Investigations  of  these 
phenomena  are  at  present  the  only  feasible  approach  toward  the  understanding 
of  the  processes  Involved. 

Modulation  of  a  wave  train  produced  by  Internal  waves  or  currents  Is  a 
strong  Interaction  effect  which  has  been  observed  In  the  laboratory,  and 
possibly  by  synthetic  aperture  radar  onboard  satellite.  More  explicitly,  a 
wave  train  or  group  can  Interact  with  a  long  wave,  and  has  a  spectral  signature 
spreading  over  a  range  of  frequencies  proportional  to  the  long  wave  speed  and 
approximately  proportional  to  the  long  wave  slope.  The  effect  of  internal 
waves  in  the  form  of  wavy  current  distribution  can  be  an  Integrated  part  of 
these  long  wave- short  wave  Interactions  without  being  smeared,  to  the  surprise 
of  almost  everyone. 

Microscale  surface  properties  have  time  spatial  densities  that  respond 
rapidly  to  energy  exchanges  with  short  gravity  waves.  Recent  theoretical  and 
laboratory  studies  indicate  that  these  densities  will  not  have  a  state  of 
saturation.  Therefore,  a  remote  sensor  estimating  the  spatial  strength  of  the 
return  signal  should  be  able  to  measure  distributions  and  variations  In  these 
densities. 

6.3  Surface  Effects  due  to  Bathymetry,  Subsurface  Mesoscale  Features, 

Internal  Waves,  Subsurface  Ekman  Layer  and  Breaking  Waves,  and  their 

Relationship  to  the  Subsurface  Processes 

As  ocean  surface  waves  propagate  from  deep  to  shallow  water  with  variable 
depth,  significant  changes  occur  In  wave  spectral  characteristics.  These 
changes  are  the  combined  effects  due  to  physical  dynamic  processes  which  may 
or  may  not  be  governed  by  the  conservation  of  energy.  The  physical  dynamic 
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processes  whose  total  wave  energy  densities  are  conserved  Include  those  processes 
such  as  wave  refraction,  wave  shoaling,  and  wave-wave  Interaction.  The  physical 
dynamic  processes  whose  total  wave  energy  densities  are  not  conserved  can  be 
those  Involved  In  energy  dissipation  and  generation.  Wave  energy  dissipation 
mechanisms  Include  percolation,  wave- Induced  bottom  motion,  wave  attenuation, 
bottom  friction,  whltecaps,  and  wave  breaking. 

The  physical  dynamic  process  Involved  In  deriving  surface  effects  due  to 
bathymetry  Is  modelled  by  the  radiative  transfer  equation  for  the  wave  spectrum. 
The  source  function  of  this  equation  Includes  those  terms  derived  for  mechanisms 
pertinent  to  the  physical  process.  Considerable  research  Is  still  required  to 
properly  model  these  source  terms. 

The  results  of  numerical  multi-layer  model  studies  Indicate  that  sub¬ 
surface  mesoscale  features  such  as  an  eddy,  a  predominant  barocllnlc  feature, 
should  exhibit  a  surface  effect  {or  signature)  in  the  form  of  deviation 
(dynamic  topography)  from  the  local  geold.  Satellite-borne  altimeter  measure¬ 
ments  have  been  abl*»  to  observe  this  deviation  Induced  by  the  Gulf  Stream,  a 
predominant  geostrophlc  barotroplc  feature,  and  also  to  confirm  the  existence 
of  deviations  due  to  eddie*.  Through  the  surface  effects  due  to  wave-current 
Interaction  (strong  interaction)  the  surface  current  boundary  conditions  of 
these  subsurface  mesoscale  features  may  further  enhance  the  possibility  of 
resolving  their  three-dimensional  density  structures  to  the  first  internal 
wave  mode  if  the  local  stability  (Brunt-Valsala)  frequency  history  is  known. 

The  Internal  waves  and  their  surface  effects,  discussed  in  Section  4.2, 
can  be  the  sources  of  the  surface  signatures  due  to  subsurface  shear  current, 
layered  density  anomaly,  bathymetry,  and  can  be  used  to  evaluate  heat  storage. 

(in  the  mixed  layer).  The  propagation  speed  of  the  internal  wave  packets, 
obtained  by  knowing  the  spacing  of  those  internal  wave  packets  or  through  the 
internal  wave-surface  wave  interactions,  depends  upon  the  density  anomaly  and 
depth  of  the  upper  mixed  layer. 

In  addition  to  surface  distortion  of  the  wave  field,  caused  by  Internal 
waves  through  strong  Interactions,  the  surface  wind-driven  Ekman  layer  will 
show  signatures  of  variations  of  Internal  vortlcity  through  horizontal 
convergence,  and  may  show  signs  of  bottom  Ekman  layer  and  free  Internal  shear 
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layer  variations  through  several  surface  signatures.  It  has  been  shown  that 
Ekman  layers  develop  rolls;  such  rolls  may  show  In  the  surface  Ekman  layer  as 
Langmuir  circulations.  However,  the  rolls  or  the  strlatlons  observed,  from 
IANDSAT  Images  and  SEASAT-1  Imagery,  seem  to  show  their  orientations  are 
parallel  to  the  velocity  field  beneath  the  Ekman  layer,  rather  than  at  an 
angle  with  the  gaostrophlc  flow.  These  strlatlons  may  be  due  to  a  periodic 
structure  In  the  Ekman  layers  at  the  bottom  or  the  Ekman  layers  at  shear 
Interfaces  beneath  the  surface.  The  direction  of  isobaths  at  the  bottom  or  of 
Isobaths  In  the  free  shear  layer  that  exist  below  the  surface  current  can  be 
deduced. 

Of  practical  importance  to  the  first  order  approximation,  a  recent 
approach  to  simulate  a  wide  variety  of  both  surface  and  subsurface  dynamic 
phenomena  has  been  developed  for  evaluating  the  Influence  of  breaking  waves  on 
the  dynamics  of  the  upper  ocean.  Based  upon  the  balance  of  energy  between 
wind  and  waves,  the  influence  was  found  to  be  dependent  of  only  one  parameter, 
the  significant  slope.  This  parameter  is  defined  as  the  ratio  of  the  root 
mean  square  wave  amplitude  to  the  dominant  wavelength  of  the  random  ocean 
sut  race  waves.  The  energy  released  from  the  breaking  waves  Is  treated  as  the 
sole  source  of  turbulence  energy.  Various  models  can  be  proposed  to  calculate 
other  parameters  such  as  wave  attenuation  rate,  equivalent  eddy  viscosity, 
wave  induced  drag  coefficient,  surface  drift,  white-capping  percentage,  and 
mixing  efficiency  for  the  mixing  layer  model. 

Their  treatment  of  using  breaking  waves  as  the  turbulence  energy  source 
alone  may  be  perceived  only  as  the  first  order  of  approximation  on  these 
parameters  or  the  lower  bound  of  their  real  values.  The  quantitative  values 
of  all  these  parameters  agree  with  the  conventional  adapted  values.  Since  the 
significant  slope  is  a  measureable  quantity  by  remote  sensors,  this  new  approach 
deserves  further  investigation. 
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